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PACS. 45. 70. -n - Granular systems. 

PACS. 81. 70. Ha - Testing in microgravity environments. 

PACS. 45.50.-j - Dynamics and kinematics of a particle and a system of particles. 



Abstract. - We report an experimental study of a dilute "gas" of inelastically colliding 
particles excited by vibrations in low gravity. We show that recording the collision frequency 
together with the impulses on a wall of the container gives access to several quantities of 
interest. We observe that the mean collision frequency does not scale linearly with the number 
N of particles in the container. This is due to the dissipative nature of the collisions and is also 
directly related to the non extensive behaviour of the kinetic energy (the granular temperature 
is not intensive). 



Introduction. - Since Maxwell and Boltzmann predictions in 1860-70, it is well-known 
that molecules of a gas move erratically with a Gaussian velocity distribution, as experimen- 
tally verified later [1]. This allows thermodynamic and transport properties of molecular gases 
to be described. However, these kinetic theory results do not hold if the particle interactions 
are dissipative or depend on their velocities (e.g., in relativistic plasmas [2]). A well-known 
example of dissipative gas is the granular gas (see [3,4] for a recent collection of papers). Since 
collisions between granular particles are inelastic, a continuous input of energy (by vibrating 
a piston or the container) is required to reach a nonequilibrium steady state. In this regime, 
granular matter sometimes seems to behave like an molecular gas in which particles follow er- 
ratic motions, but several experiments have displayed striking different properties: Instability 
of the homogeneous density state leading to cluster formation [5-7], non-Gaussian nature of 
the velocity distribution [8], anomalous scaling of the pressure [7,9]. These effects have been 
also numerically simulated and some of them have been theoretically understood [10]. 
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In this paper, we report a 3D experiment of a dilute granular medium fluidized by sinusoidal 
vibrations in a low gravity environment. The motivation for low gravity is to achieve an 
experimental situation in which inelastic collisions are the only interaction mechanism, and 
where only one "input" variable (the inverse vibration frequency) has the dimension of time [6] . 
This eliminates possible resonances between the time of flight of a particle under gravity 
and the period of vibration. The aim is to observe new phenomena which result from the 
inelasticity of the collisions, thus absent in molecular gases. We first study the scaling of the 
collision frequency with a container wall with respect to the vibration velocity, V, and the 
particle number, N. We also measure the time lag distribution between successive collisions 
with the wall and the impulse distribution. We show in particular that two measurements 
display significant differences from the behaviours observed in molecular gases: the scalings 
of the collision frequency and of the particle impulse distribution with TV. The scaling of 
the granular temperature with V has been extensively investigated [3], but there exists only 
one 2D experiment for the scaling with N [11]. We emphasize that velocity distributions in 
granular gases have been measured so far only for nearly 2D geometries. In the 3D case, it is 
much easier to measure the distribution of impact velocities at a boundary as done here. This 
measurement involves a similar information content and can be easily compared to molecular 
dynamics simulations [12]. 
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Fig. 1 - Typical time recording of the force sensor [impulse response I(t)] during 10 periods of vibration 
showing 106 collisions. Inset: zoom of this signal during 1.3 ms showing two detected collision peaks 
(o) and the typical damping time of the oscillatory response of the sensor. The parameters of vibration 
are: N = 12, / = 40 Hz, A = 1.96 mm (not listed in TableHJl. 

Experimental setup. - A fixed transparent Lexan tube, D = 12.7 mm in inner diameter 
and L = 10 mm in height, is filled with N steel spheres, d = 2 mm in diameter. Experiments 
have been performed with N = 12, 24, 36 and 48 respectively corresponding to n = 0.3, 0.6, 
0.9 and 1.2 particle layers at rest (packing fraction from 0.04 to 0.18). A piezoelectric force 
sensor (PCB 200B02), 12.7 mm in diameter, is fixed at the top of the cell in order to record 
the particle collisions with the upper wall. A piston made of duralumin, 12 mm in diameter, 
is driven sinusoidally at the bottom of the cell by an electromagnetic shaker. The frequency / 
is in the range 40 to 91 Hz and the maximal displacement amplitude A is varied from 0.4 to 2 
mm. The vibration parameters during the time line are listed in TableHJ Vibration amplitudes 
are measured by piezoelectric accelerometers (PCB 356A08) screwed in the shaft in a triaxial 
way. Typical output sensitivities in the vibration direction and in the perpendicular directions 
are, respectively, 0.1 and 1 V/g, where g = 9.81 m/s 2 is the acceleration of gravity. Typical 
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force sensor characteristics are a 11.4 mV/N output sensitivity, a 70 kHz resonant frequency, 
and a 10 /xs rise time. Low gravity environment (about ±5 x 10 g) is repetitively achieved by 
flying with the specially modified Airbus A300 Zero-G aircraft through a series of parabolic 
trajectories which result in low gravity periods, each of 20 s. An absolute acceleration sensor 
allows the detection of the low gravity phases and the automatic increment of the vibration 
parameters after each parabola. The output signals of force, respectively accelerations, are 
stored on a computer on 16 bits at a 2 MHz sampling rate, respectively, on 12 bits at 10 kHz. 
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Table I - Vibration parameters during each parabola. V = lirAf and Y = A-K 2 Af 2 /g are respectively 
the maximal piston velocity and the dimensionless acceleration. The number of collisions N w on the 
sensor is detected during 8 = 16 s of low gravity to avoid transient states. 



Detection of collisions. - A typical time recording of the force sensor shows a succession 
of peaks corresponding to particle collisions, as displayed in Fig.Q]for 10 periods of vibration. 
Bursts of peaks roughly occur in phase with the vibration but the number of peaks in each 
burst and their amplitude are random (see Fig. . A peak corresponds to the collision of 
a single sphere, which leads to an almost constant impact duration from 5 to 6 fis for our 
range of particle velocities v (assumed of the order of V). Indeed, the Hertz's law of contact 
between a sphere of radius R and a plane made of same material, leads to a duration of 
collision r = YR/v 1 ^ 5 , where Y = 6.9 x 10~ 3 (s/m) 4 / 5 for steel [13]. The signal recorded 
by the sensor corresponds to an impulse response, I(t). Each peak due to a collision is thus 
followed by an oscillatory tail at the sensor resonance frequency (roughly 100 kHz) damped 
over 500 //s (see inset of Fig.QJ. A thresholding technique is applied to the signal in order to 
detect the collisions. We have to discard a time interval of 100 /is around each detected peak 
in order to avoid counting the first maxima of each oscillatory tail as additional collisions. 
Thus an additional weak collision occuring in the oscillatory tail due to the previous one may 
be missed by our detection process. However, the discarded time interval is small compared 
to the mean time lag (a few ms) between successive collisions if their statistics is assumed 
Poissonian (see below). Consequently, the probability of possibly discarded collisions is small. 

Collision frequency scaling. - The number of collisions N w with the top wall (i.e., the 
sensor) is obtained by the previous thresholding technique, for each parameter listed in Table 
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HI during 9 = 16 s of low gravity to avoid possible transient states. For a fixed number of 
particles, N, Fig. shows that the collision frequency, v w = N w /9, is proportional to the 
maximal piston velocity, V, for 0.1 < V < 0.5 m/s. As also shown in Fig. [2} v w oc VN a , 
with a = 0.6 ± 0.1 for our range of TV. This result strongly differs from the kinetic theory of 
molecular gases for which v w varies linearly with N. It cannot be explained either in the very 
dilute limit (Knudsen regime). Indeed, assuming that each particle mostly collides with the 
boundaries of the container and does not interact with others, leads to v w oc VN/[2(L — d)}. 
Therefore particles do interact significantly with each other through inelastic collisions. We 
will show below that this anomalous scaling is a consequence of the dissipative nature of 
collisions. Note also that this scaling law with N has been recently recovered in 2D numerical 
simulations [12]. It thus appears to be a robust and generic behaviour of granular gases as we 
will explain below. 
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Fig. 2 - Collision frequency rescaled by the number of particles, i/ w /N ' 6 , as a function of V for 
N —12 (□) and (■); 24(0); 36(y); 48(o). ■-marks are from a previous set of experiments at fixed 
N = 12 for 15 different velocities which are not listed in Table U Solid line corresponds to the fit 
Vw/N - 6 = V/lo where l ^ 5.9 mm. 

Time lag distribution. - The probability density functions (PDF) of the time lag At 
between two successive collisions with the top wall is displayed in Fig. for 4 different values 
of N and various parameters of vibration. These PDFs are found to decrease exponentially 
with At and to scale like V for our range of V. This exponential distribution for the time 
lag statistics is the expected one for Poissonian statistics. As already shown for the data of 
Fig. |21 these PDFs can be collapsed by the TV - 6 rescaling. We also observe in Fig. that 
even the largest values of AtV are smaller than L. In our range of N, the Knudsen number, 
K = l/L, is in the range 0.1 — 1, where I is the mean free path, I = £l/(Nnd 2 ). We are thus in 
a transition regime from a Knudsen regime to a kinetic regime. It corresponds to a crossover 
between the very dilute regime for which each particle mostly collides with the boundaries (I 
of order L independent of N) , to the kinetic regime (I inversely proportional to N) . Finally, if 
the amplitude of vibration, A, is not negligible with respect to L (i.e., A/L > 0.17), the time 
lag distributions are no longer exponential (not shown here). 

Impulse distribution. - The PDFs of the maxima I of the impacts recorded in Fig. 
are displayed in Fig. for 4 different values of N and various parameters of vibration. Note 
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Fig. 3 - Probability density functions of the time lag At between two successive collisions rescaled 
by V (V = 2ivAf), for N = 12, 24, 36 and 48 particles, and for different vibration parameters (see 
symbols in Table QJ. 




Fig. 4 - Probability density functions of the impulse I of the impacts on the sensor rescaled by V, 
for different vibration parameters. Symbols are the same as Fig- El 
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that the low impulse events are not resolved because of noise. Indeed, we expect that the 
PDFs vanish for 1 = 0. We first observe that they scale like V, for our range of V (see Fig. 
QJ. Second, they display exponential tails with a slope increasing with the particle number 
TV. Third, the PDF for different values of TV can be roughly collapsed when / is scaled like 
V/N 13 with ^0.8± 0.2. This shows that the mean particle velocity v near the wall scales 
like, v cx V/N 13 , which gives for the granular temperature near the wall, T w cx V 2 /N 2 @ (see 
below). Finally, we observe that the shape of the distributions of / also differ from kinetic 
theory of molecular gas. They display an exponential tail instead of the Gaussian one. Note 
however, that it has been observed many times that a universal shape of the bulk velocity 
distribution do not exist for granular gases. In particular, a significant effect of side walls has 
been reported [8]. 

Discussion and concluding remarks. - The experimental results obtained on the collision 
frequency and the distribution of impulse at the wall are related and can be used to extract 
information on various quantities of interest. Indeed, keeping only quantities that depends 
on V and TV, we have v w cx p w v cx p w I where p w is the particle density at the wall. Thus, 
we get p w cx N a+I3 and T w cx V 2 N~ 213 . The density of particles close to the wall opposite to 
the piston increases faster than TV (a + (3 = 1.4 ± 0.3) because the density gradient becomes 
larger when TV increases. Indeed, it is well known that the granular temperature decreases 
away from the piston because of inelastic collisions. Thus, the density has to increase in order 
to keep the pressure P constant (in zero gravity environment). We also have from the state 
equation for a dilute gas, P cx p w T w cx V 2 N a ~ 13 . 

The dependence on TV of the collision frequency can be understood as follows: we have 
v w cx p w v. In a molecular g ./V and v is fixed by the thermostat and does not depend 

on TV, thus we have v w cx TV. In a granular gas, v or the total kinetic energy E are determined 
from the balance between the injected power by the vibrating piston and the dissipated one 
by inelastic collisions. For a dilute gas with a restitution coefficient r very close to 1 such that 
the density is roughly homogeneous (p w cx TV), E does not depend on TV, and thus v cx l/VTV. 
Thus, we get v w cx \/N [12]. In the present experiment, both the scaling of the density and 
the one of the mean velocity differ from this limit case, because gradients of density and 
granular temperature cannot be neglected. However, the prediction in the limit r ~ 1 gives 
a good approximation to the observed scaling of the collision frequency. This shows that a 
granular gas driven by a vibrating piston strongly differs from a molecular gas in contact with 
a thermostat even in the limit r ~ 1. 
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